Unconventional London penetration depth in Ba(Feo.93Coo.o7)2As2 single crystals 
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The London penetration depth, A(T), has been measured in several single crystals ol 
Ba(Feo.93Coo.o7)2As2. Thermodynamic, electromagnetic, and structural characterization measure- 
ments confirm that these crystals are of excellent quality. The observed low temperature variation 
of A(r) follows a power-law, AA(r) ~ T" with n — 2.4 ± 0.1, indicating the existence of normal 
quasiparticles down to at least 0.02Tc. This is in contrast to recent penetration depth measure- 
ments on single crystals of NdFeAsOi-i^F^ and SmFeAsOi-a;Fa;, which indicate an anisotropic but 
nodeless gap. We propose that a more three-dimensional character in the electronic structure of 
Ba(Feo.93Coo.o7)2As2 may lead to an anisotropic s— wave gap with point nodes that would explain 
the observed A(r). 

PACS numbers: 74.25.Nf,74.20.Rp,74.20.Mn 



The discovery of 
LaFeAs(Oi_^F^), with 



superconductivity in 
23 K 1], and above 50 



K in SmFeAs(0i_3:F2,) 0] has lead to a flurry of exper- 
imental and theoretical activities aimed to understand 
the fundamental physics governing this new family of 
superconductors. Much effort has been put forth to 
characterize the oxypnictide compounds RFeAsOi-a^Fj, 
(R=rare earth, "1111" in the following), as well as 
oxygen-free A(Fei_a;T2:)2As2 (A — alkaline earth and 
T = transition metal, "122" in the following). The 
major questions yet to be answered are the symmetry 
of the superconducting gap and its universality among 
all Fe-based pnictides. High transition temperatures, 
high upper critical fields, and the existence of a unique 
structural element, Fe-As layers, prompt a comparison 
with the high-Tj, cuprates. 

It has now been well established that the cuprates ex- 
hibit d-wave symmetry of the order parameter In 
Fe-based pnictides the debate is wide open. Using dif- 
ferent techniques, several groups have arrived at differ- 
ent conclusions, such as gapped versus nodal, or sin- 
gle versus multigap superconducting states. Tunneling 
measurements showed an unconventional order param- 
eter with nodes [1], two superconducting gaps and a 
pseuodogap as well as conventional s-wave BCS be- 
havior [6| . Measurements of specific heat have also shown 
both conventional [tI and unconventional [11 behavior. 
Angle-resolved photoemission (ARPES) data are still in 
disagreement regarding the gap anisotropy, amplitude, 
and coinparison to theoretical electronic structure calcu- 
lations 

As far as penetration depth is concerned, measure- 
ments reported for the 1111 system are consistent with 
regard to the overall gap topology. Microwave cavity per- 
turbation in PrFeAsOi_:r [111 and muon spin relaxation 
(^SR) in both LaFeAsOi-^r ,12] and SmFeAsOi_:,F^ ^ 
are consistent with a nodeless order parameter. Single 
crystal measurements of the London penetration depth 



using a tunnel diode resonator (TDR) in SmFeAsOi-ajF^, 
[l^ and NdFeAsOi_a;Fa; [l^ have both found exponen- 
tial behavior in A(T 0). The superfluid density has 
been successfully fit to a single moderately anisotropic 
gap [ll] and a two gap model without nodes [l3]- How- 
ever, the question of the order parameter symmetry in 
the 122 iron pnictides is still open. Over the past two 
decades much attention has been devoted to the cuprates, 
where the band structure is essentially two-dimensional 
and line nodes in the gap result in a A(T) ~ T behavior 
in the clean limit [sl , changing to ~ by the impurities 
[l^ . Still, other materials may have nodes in a three- 
dimensional gap. There is evidence f or p oint nodes in 
UBei3 UPtg [3, PrOs^bi2 [ll,|20|, and possibly 



nonmagnetic borocarbides |2l| . 

In this Letter, we focus on the Ba(Feo.93Coo.o7)2As2 
member of the 122 family, for which high quality single 
crystals are available 2^, 23, The behavior in the 
mixed state is similar to the high Tc cuprates, such as Y- 
Ba-Cu-0 and Nd-Ce-Cu-0 [2J], and one could expect to 
find unconventional behavior in other properties as well. 

The parent compound BaFe2As2 has been studied in 
detail elsewhere [25| . Single crystals of superconducting 
Ba(Fei_2;Co2:)2As2 were grown out of FeAs flux using 
standard high temperature solution growth techniques. 
The FeAs and CoAs powders were mixed with Ba accord- 
ing to the ratio Ba:FeAs:CoAs=l:3.6:0.4, placed into an 
alumina crucible, and then sealed in a quartz tube un- 
der an approximately 1/3 atm of argon gas. The sealed 
quartz tube was slowly heated to 1180 °C, held for two 
hours, and then cooled to 1000 °C over the course of 30 
hours. Upon reaching 1000 °C the FeAs was decanted 
from the single crystals. The crystal dimensions can be 
as large as 12 x 8 x 1 mm'^. 

Crystals selected from several different batches have 
been extensively characterized by various techniques. 
Elemental analysis has shown that the actual Co 
concentration (averaged over six measurements) is 
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FIG. 1: (Color online) Characterization of 

Ba(Feo.93Coo.o7)2As2 single crystals. (a) Longitudinal 
x-ray scans through the position of the (008) reflection for 
indicated temperatures, (b) similar (^ ^ 0) scans through 
the position of the (118) reflection, (c) In-plane resistivity 
(open circles), dc magnetization (filled circles) and in- plane 
penetration depth (squares). 



Co/(Co+Fe)=7.4 % ±0.2 %. Powder x-ray diffraction on 
ground crystals has yielded tetragonal lattice constants 
of a = 3.9609 ± 0.0008 A and c = 12.9763 ± 0.004 A. 
Single crystal x-ray diffraction measurements have found 
no evidence of a tetragonal to orthorhombic structural 
transition, as shown in Fig. [Ua)-(b), and the compara- 
tive analysis of the peak widths indicates homogeneous 
doping. Figure[T](c) shows in-plane resistivity, dc magne- 
tization, and TDR penetration depth. Tc, as defined by 
zero resistivity, coincides with the onset of magnetization 
and is close to the shoe of the TDR data. 

Measurements of the London penetration depth were 
performed using a tunnel diode resonator technique 

m 



271 l28j . A properly biased tunnel diode acts as 



an ac power source for the LC circuit that promotes 
spontaneous oscillations at a natural frequency, 27r/o = 
where Lq is the inductance of the coil with- 
out a sample inside. The sample, mounted on a sap- 
phire rod, is inserted into the coil without touching it. 
Throughout the measurement, the circuit is stabilized 
at 5.0 K ± 3 mK to ensure a stability in the resonance 
of less than 0.05 Hz over a period of several hours. In 
our setup /o ~ 14 MHz, which implies that the fre- 
quency resolution is better than 5 parts in 10^. This 
translates into ^ 1 A resolution in AA for mm-sized 
crystals. For samples much smaller than the coil (small 
filling factor), the data are noisier, as can be seen in 
the data for NdFeAsOi_i,F^ in Fig. [H Due to the dia- 



magnetic susceptibility x of the sample, the inductance 
changes to a new value, L, and the resonant frequency 
shifts accordingly. A/ = / — /o ~ —GAt:x{T), where 
G = foVs/2Vc (1 — N) is the calibration constant, Vs is 
the sample volume, Vc is the effective coil volume and N 
is the demagnetization factor. For a superconductor in 
the Meissner state, Anx = A/i?tanhi?/A — 1, where the 
effective dimension R takes into account the actual sam- 
ple shape [i^]. While in the coil, the sample experiences 
a ^ 10 mOe ac excitation field, much smaller than the 
first critical field, Hd, and therefore the London pene- 
tration depth is measured. The calibration constant G 



has been determined by two different techniques [27|, |28( . 
First, the sample was physically pulled from the coil at a 
low temperature to determine the full frequency change 
due to its presence. Second, the normal state skin depth 
was calculated from the measured resistivity just above 
Tc, Fig. [TJc). Both approaches yield the same result, 
which also rules out inhomogeneity of our samples. The 
TDR technique precisely measures changes in the pene- 
tration depth, but its absolute value is more difficult to 



obtain [28|. To estimate A(0), we determined the lower 



critical field, Hd, by measuring the M{H) loops while 
increasing the maximum field until nonlinearity due to 
vortices appeared. Using Hd = $o/(47rA'^)(ln A/^ + 0.5), 
with a value for the coherence length of ^ = 2 nm, as 
estimated from the upper critical field, we have obtained 
A(0) « 208 nm. This value is similar to other reports, 
where A(0}_R:! 254 nm for La-1111 ^ and 190 nm for 
Sm-1111 [13| from /xSR measurements. 

The low temperature behavior of AA(T) = \{T) — 
A(0) and the superfluid density, ps{T) = (A(0)/A(^))^ 
are commonly used to determine the symmetry of 
the superconducting pairing state [29| . In the case 
of a fully gapped Fermi surface, AA(T)/A(0) ~ 
^ttA (0) /2kBTexp (-A (0) /ksT), which is also true for 
the cases of an anisotropic gap and two-gap supercon- 
ductivity if one allows A(0)/rc to be a free parame- 
ter. For the isotropic s-wave BCS case, this form is 
applicable for T < Tc/3. This is consistent with the 
data taken in our TDR system for the conventional BCS 
superconductor Nb, Fig. [2l An wave behavior, al- 
beit with an anisotropic gap with smaller amplitude, is 
found in NdFeAsOi-^^F:, [ll]. Fig. H In the case of 
d-wave pairing, AA(T)/A(0) « r[21n2/aA(0)], where 
a = A-i (0) \dA (0) [s^. This linear behav- 

ior is observed in our measurements on clean BSCCO- 
2212 crystals, Fig. [21 The above equations are appli- 
cable only for clean cases. Impurity scattering does 
not affect the isotropic s-wave gap but it does sup- 
press the d-wave gap. The superfluid density, however, 
does change in both cases, as shown in FigO In the 
case of d-wave pairing, scattering leads to AA(T) ~ 
[l^ . For an s-wave superconductor in the dirty limit, 
Ps{T) = A(r)/A(0)tanh [A(T)/2fcBT]P. 

The observed low-temperature behavior of A(T) in 




FIG. 2: (Color online) Comparison of the low temperature 
AA(T) measured in Ba(Feo.93Coo.o7)2As2 with the known d- 
wave (BSCCO) and s-wave (Nb) superconductors as well as 
NdFeAsOi-^F:!,. Inset: AA vs. {T/Tcf for three differ- 
ent samples of Ba(Feo.93Coo.o7)2As2 from different batches. 
(Curves are offset for clarity). 
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FIG. 3: (Color online) Superfluid density, ps, as a func- 
tion of normalized temperature for Ba(Feo.93Coo.o7)2As2 us- 
ing A(0) = 200 nm and 300 nm. Also shown for comparison 
are calculated curves for clean and dirty s-wave and d-wave 
pairings as well as best fit of NdFeAsOi-^Fa; compound (cir- 
cles) from ReflTsl. 



Ba(Feo.93Coo.o7)2As2 is clearly non-exponential, Fig. 
The inset of Fig. [2] shows the data plotted as a function 
of (T/Tc)^ for samples from different batches all having 
different sizes and aspect ratios. The middle curve corre- 
sponds to the smallest sample and therefore the small- 
est filling factor. The best fit to A(T) - T" yields 
n = 2.4±0.1, where the error reflects the fitting of differ- 
ent curves. For comparison, the penetration depth for the 
NdFeAsOi-^F^: crystal, plotted as open circles in Fig. [21 
is fiat below about 0.2Tb, providing evidence for a fully 
gapped Fermi surface [15|. 

To gain a better understanding of the pairing state, in 
Fig. [3] we plot Ps{T) along with the known s— and d— 
wave behavior and the best fit for NdFeAsOi-^^F^, [3]. 
The data for the Ba(Feo.93Coo.o7)2As2 samples are plot- 
ted for A(0) = 200 nm and 300 nm. The inset of Fig. [3] 
zooms into the low temperature region showing that 
the observed non-exponential behavior persists down to 
0.02Tc. As described above, the value of A(0) « 200 
nm was found from Hd and is consistent with the liter- 
ature [H, The curve for A(0) = 300 nm produces 
Ps{T) closer to standard curves, but no A(0) can change 
the non-exponential low temperature behavior. At in- 
termediate temperatures, the experimental Ps{T) shows 
significant departure from the d-wave and s-wave curves. 
Fig. EI Possible reasons for this behavior could be signif- 
icant gap anisotropy or multi-gap pairing with different 
gap amplitudes. We note that for NdFeAsOi-j^F^, crys- 
tals, a similar behavior in the intermediate temperature 
region has been found (Fig. [3] open circles) and explained 
in terms of an anisotropic s-wave gap [l5|. However, this 
does not lead to the low-temperature behavior. 



A possible explanation for the unusual exponent in 
A ^ T" with n > 2, is the existence of a three- 
dimensional (3D) gap with point nodes. In such a case, 
a quadratic behavior of A(T) is predicted and observed 
in the clean case [13; 113 • With impurities, the exponent 
n increases, consistent with the observed n = 2.4. Sim- 
ilarly, in the 2D d-wave case T— linear behavior of the 
penetration depth changes to [3]. Since a three- 
dimensional portion of the Fermi surface is a crucial 
ingredient for the point-node scenario, electronic struc- 
ture calculations. Fig. 21 were performed within the full- 
potential linearized plane wave (FLAPW) approach [3^ 
using the local density approximation (LDA) [33]. Both 
relaxed and experimental As positions were used with the 
experimental lattice constants for the BaFe2As2 [1^ and 
Ba(Feo.g3Coo.o7)2As2 (see above). A virtual crystal ap- 
proximation was employed with Fe replaced by a virtual 
atom having Z=26.07 to simulate 7.4% Co doping. The 
Fermi velocities, calculated using the BolzTrap [3^] pack- 
age, are very sensitive to the As positions, a situation 
similar to calculations of magnetic properties [s^, 37 1. 



The calculated coordinate, Zyis=0.341, is comparable to 
previous calculations, z^s=0.342 [37], but is significantly 
lower than the experimental value of z^s=0.355. This 
downshift by 0.16 A increases the band dispersion along 
the z direction. The calculations with relaxed zas give 
7^ = (^k) / {''^z) — 3 for pure compound and 7^ = 2.7 
for Ba(Feo.93Coo.o7)2As2. With the experimental zas, 
the difference is much larger, 7^ = 12.1 and 9.0, for pure 
and doped materials, respectively. 

The fact that the superconducting gap functions in 
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FIG. 4: (Color online) (a) Fermi surface of non- 
superconducting BaFe2As2 and (b) Ba(Feo.93Coo.o7)2As2 
with experimental As positions. Note significant enhance- 
ment of a 3D character around F, Z and X. 

NdFeAs(Oi_a;Fj;) and Ba(Feo.93Coo.o7)2As2 behave so 
differently may imply that the two systems are gov- 
erned by truly distinct pairing mechanisms. However, 
both systems share the Fe-As planes as key structural 
building blocks. Electronic structure calculations demon- 
strate that the states at the Fermi energy are domi- 
nated by Fe-3c? states. Thus, we discuss the implica- 
tions of our findings under the premise that there is a 
common pairing mechanism and a common pairing sym- 
metry in both materials. Our earlier experimental re- 
sults for NdFeAsOi_2:F2: rule out a d-wave pairing state. 
Previously, we argued that a likely interpretation of the 
moderate anisotropy of the superconducting gap, A ((^), 
[^, can naturally be understood as a result of an elec- 
tronic interband pairing interaction fsO^ that changes the 
sign of the gap. In this scenario, anisotropic s-wave pair- 
ing occurs with the nodes located between Fermi surface 
sheets [i^l . Then, A (ip) ^ for all if and is slightly 
smaller along the diagonals of the Brillouin zone con- 
necting F and X points. As shown in[4l the Fermi sur- 
face sheets of Ba(Feo.93Coo.o7)2As2 are more 3D than 
those of the parent BaFe2As2, and even more so than 
NdFeAs(0i_j;F3;) [11]. The "bellies" of the sheets around 
the X— points in the Brillouin zone corners reach out to- 
ward the nodes and a vanishing gap can easily occur at 
a point or small ellipse close to the "belly button" . As in 
NdFeAs(0i_2;F2;), this reduction of the gap is expected 
to occur along the directions connecting the F and X 
points. Such an anisotropic s-wave state with nodes im- 
plies that the existence of nodes depends sensitively on 
the Fermi surface shape, in sharp contrast to d- or p-wave 
states where the nodes are enforced by the symmetry. 
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